APPROXIMATION OF FOURIER INTEGRAL OPERATORS BY 

GABOR MULTIPLIERS 

ELENA CORDERO, KARLHEINZ GROCHENIG AND FABIO NICOLA 

Abstract. A general principle says that the matrix of a Fourier integral oper- 
ator with respect to wave packets is concentrated near the curve of propagation. 
We prove a precise version of this principle for Fourier integral operators with a 
w3 ' smooth phase and a symbol in the Sjostrand class and use Gabor frames as wave 

packets. The almost diagonalization of such Fourier integral operators suggests 
a specific approximation by (a sum of) elementary operators, namely modified 
Gabor multipliers. We derive error estimates for such approximations. The 
methods are taken from time-frequency analysis. 
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1. Introduction 



< 

r-| ! A fundamental principle expressed by Cordoba and Fefferman [10] says that 

c^ I Fourier integral operators map wave packets to wave packets. Furthermore each 

wave packet is transported according to the canonical flow in phase space that is 
associated to the operator. 

Rigorous versions of this principle were proved for various types of Fourier inte- 
^ I gral operators and wave packets. Cordoba and Fefferman used generalized Gaus- 

sians as wave packets and considered Fourier integral operators with constant sym- 
O ■ bol. In the last decade the transport of wave packets by Fourier integral operators 

^ ! was investigated for curvelet-like frames by Smith ^29j, for curvelets by Candes and 

^ I Demanet [3], for shearlets by Guo and Labate [15], and more recently for Gabor 

frames by Rodino and two us [H]- 

In this paper we investigate Fourier integral operators with a "tame" phase and 
non-smooth symbols with respect to time-frequency shifts and Gabor frames. Fol- 
lowing [1] , we will assume that the phase function $ on R^*^ satisfles the conditions 
(i) $ e C°°(M2'^), (ii) 9^$ G L°°(M2'^) for a > 2, and (iii) mi^^^^^2, \ det 9^ ,^$(x, r])\ > 
S > 0. For brevity, we call $ a tame phase function. Then the Fourier integral 
operator T with phase function $ and the symbol a is formally deflned to be 

(1) Tf{x)= [ e'^''"^''''^a{x,v)mdv- 



Fourier integral operators with a tame phase were investigated by many authors, 
both in hard analysis [2] and in time-frequency analysis [5], [H [9]. They arise the 
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study of Schrodinger operators, for instance, for the description of the resolvent 
of the Cauchy problem for the Schrodinger equation with a quadratic type Hamil- 
tonian [H [221 123]. In our preceeding work [SI El [7] we have adopted the point of 
view that Fourier integral operators with a tame phase are best studied with those 
time-frequency methods that correspond to a constant geometry of the wavepack- 
ets. Precisely, let z = {x, r]) G M^'^ be a point in phase space and 

7r{z)f{t) = e'-''^-'f{t-x), teR" 

be the corresponding phase space shift (or time-frequency shift). Roughly speaking, 
these wave packets correspond to a uniform partition of the phase space. 

In this paper we continue the investigation of the matrix of a Fourier integral 
operator with respect to (frames of) phase space shifts. A first version of a time- 
frequency analysis of Fourier integral operators was proved in [8J. The formulation 
is in the spirit of Candes and Demanet and is intriguing in its simplicity. 

Theorem 1.1. Assume that ^ is a tame phase function and x is the corresponding 
canonical transformation ofM.'^'^. If g G iS(]R'^) and if a belongs to the Sobolev space 
W^{M.^'^), then there exists a constant Cn > such that 

(2) \{T7c{fi)g,n{X)g)\<CN{x{f^)-X)-''', VA,/i eR''. 

Our first goal is to sharpen Theorem 11.11 in several respects. 

(a) We keep the assumptions on the phase function $ (tame phase), but allow 
non-smooth symbols taken in a generalized Sjostrand class. 

(b) We enlarge the class of "basis" functions or "windows" g and extend Theo- 
rem [TTT] to hold for windows in certain modulation spaces M^{R'^). 

The appearance of modulation spaces is not longer a surprise. It is well estab- 
lished that modulation spaces are the appropriate function spaces for phase space 
analysis and time-frequency analysis with constant geometry. In particular the 
Sjostrand class and its generalizations are a special case of the modulation spaces 
and they arise as suitable symbol classes for pseudodifferential operators. The use 
of the Sjostrand class for Fourier integral operators goes back to Boulkhemair [2] 
and features prominently in recent work on pseudodifferential operators and Fourier 
integral operators (see [T^] for a survey). Furthermore, the theory of time- frequency 
analysis and modulation spaces is now well developed and offers new tools for the 
investigation of Fourier integral operators, in addition to the classical methods. 

The arguments and proofs of the extension of Theorem 11.11 are in the spirit of 
[H [18], but require a substantial amount of technicalities. A main technical result 
is a new characterization of the Sjostrand class in Proposition 13.101 This charac- 
terization is perfectly adapted to the investigation of Fourier integral operators and 
measures modulation space norms with a variable family of windows rather than a 
fixed window. 

Our second goal is to derive a discretized version of Theorem 11.11 We will show 
that a Fourier integral operator can be approximated by (sums of) elementary 
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operators of the form 

(3) MJ = M^'^'^f = J2a,{f,n{X)g)n{x'{X))9, 

aga 

where A is a lattice in M^*^, [ax] a symbol sequence on A, and x' : A — )■ A is a discrete 
version of the canonical transformation associated to a phase function $. Such 
operators are suggested by a recent approximation theory for pseudodifferential 
operators through Gabor multipliers |20]. The modified Gabor multipliers in ([3]) 
are adapted precisely to the canonical transformation. We will approximate a 
Fourier integral operator with tame phase by modified Gabor multipliers and prove 
an approximation theorem for Fourier integral operators (Theorem 15 ■4p . 

The approximation theory also yields an alternative proof for the boundedness 
of Fourier integral operators on modulation spaces. The boundedness and Schat- 
ten class properties of Fourier integral operators were studied under various as- 
sumptions on the phase and the symbol. We refer to the book [30] and the arti- 
cles [21 [TOl O [TJ |27] for a sample of contributions. 

Gabor multipliers have many applications in signal processing and acoustics [HI 
|3T] and are especially useful for the numerical reahzation of pseudodifferential oper- 
ators. For the study of Fourier integral operators we introduce a new type of Gabor 
multipliers. We hope that the modified Gabor multipliers will also prove useful for 
the numerical realization and approximation of Fourier integral operators. 

Our paper is organized as follows: In Section 2 we collect some concepts and the 
necessary facts from time-frequency analysis. Section 3 is devoted to the almost 
diagonalization of Fourier integral operators with a tame phase. We will prove 
a substantial generalization of Theorem 11.11 In Section 4 we define and study 
the modified Gabor multipliers mentioned in ([3]) and prove their boundedness on 
modulation spaces. In Section 5 we study the representation and approximation 
of Fourier integral operators by sums of modified Gabor multipliers and derive a 
quantitative approximation theorem. Further remarks are contained in Section 6. 

2. Notation and preliminary results 

We recall some notation and tools from time-frequency analysis. For an exposi- 
tion and details we refer to the book [T7] . 

(i) First we define the translation and modulation operators 

TJit) = fit-x) and M,/(t) = e^-^^VW, 
for A = {x,r]) G M^'^, / G iS'(M'^). Their composition is the time-frequency shift 

7r(A) = M^T,. 

(ii) Now consider a distribution / G iS'(M'^) and a Schwartz function g G iS(M'^)\{0} 
(which will be called window). The short time Fourier transform of / with respect 
to g is given by 



VJ{x, T]) = (/, M,T,g) = (/, 7r{X)g) = / f{t)g{t - x)e'^^''-^ dt, for x,v e 
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The last integral makes sense for dual pairs of function spaces, e.g., for f,g E 
L^(]R'^). The bracket (■, ■) makes sense for dual pairs of function or distribution 
spaces, in particular for / G iS'(]R'^) and g G iS(]R'^). 

(iii) In the sequel, we set v{x, rf) = Vs{x,r]) = {{x,r])Y = {l + \x\^ + \r]\^y/'^ and we 
denote by A4.v{^'^) the space of f-moderate weights on M^"'; these are measurable 
functions rra > satisfying m{z + Q < Cv{z)m{() for every z,( ^ ^'^■ 

(iv) Then, for 1 < p < oo and m G A^t,(M^'^) we denote by M^(R°') the space of 
distributions / G 5'(]R'^) such that 

ll/llAf^(lRd) := l|K;/||Lf„(i;2d) < oo. 

This definition is meaningful and does not depend on the choice of the window 
g G S(R'^),g 7^ 0. If fact, different g G M^iR'^) yield equivalent norms on M^ (M^) 
whenever m G Mv(R'^'^)- 

(v) Modulation spaces possess a discrete description as well. Consider a lattice 
of the form A = AZ'^'^ with A G GL{2d, R). The collection of time-frequency shifts 
Q{g,A) = {7r(A)(? : A G A} for a non-zero g G L^(M'^) is called a Gabor system. 
The set G{g, A) is a Gabor frame, if there exist constants A,B > such that 

(4) AWfWl < 5^ |(/,7r(A)^)r < 5||/||^ V/ G L'{R'). 

AeA 

Gabor frames give the following characterization of the Schwartz space iS(]R'^) and 
of the modulation spaces Mp^(R'^). Fix g G S(R'^),gj^ 0, then 

(5) / G SiR') ^ sup(l + |A|)^(/, n{\)g) < oo ViV G N, 

AeA 

(6) / G Mi(R') ^ (5^|(/,vr(A)^)|MAr)'^' < oo, 

AeA 

and the latter expression gives an equivalent norm for M^{R'^). 

Finally we recall that, ii A = B = 1 in (j4]), then Q{g,A) is called a Parseval 
frame, and (jl]) implies the expansion 

(7) / = 5^(/,7r(A)^)7r(A)(7 yfeL^R^) 

AeA 

with unconditional convergence in L'^{R'^). If / G M^(R'^) and g G M^{R'^) (m G 
A^^(M^'^)), then this expansion converges unconditionally in M^(R'^) for 1 < p < 
oo. 

(vi) Amalgam spaces: Let m G A^^(M^'^). The amalgam space W{Lp){R'^) con- 
sists of all essentially bounded functions such that the norm 

\\f\\w{L?^) = Yl ^^P \f{x + k)\m{k) <oo 

is finite. If / G W{L^) is also continuous and A is a lattice in R'^, then the 
sequence (/(A))AeA is in the sequence space i^i-^)- We will use several times the 
usual convolution relations for amalgam spaces 

(8) LI * WiLl) C W{Ll,) , 
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which hold for 1 < p < oo and all m E Aiv(M?'^). 

3. Almost diagonalization of FIOs 

For a given function / on Mf^, the Fourier integral operator (FIO in short) with 
symbol a and phase $ can be formally defined by 

(9) T/(x)= / e2-*(^''')a(x,r/)/(r7)dr7. 



To avoid technicalities we take / G iS(M'^) or, more generally, / G M^{M.'^). If 
a G L°°(M^'^) and if the phase $ is real, the integral converges absolutely and 
defines a function in L°°(M'^). 

We will consider a class of Fourier integral operators which arises in the study 
of Schrodinger type equations. This class of operators was already studied in the 
spirit of time-frequency analysis in [8]. 

Definition 3.1. A phase function ^{x,7]) is called tame, z/$ satisfies the following 
properties: 

(i) $ G C°°(M2'^); 
(ii) for z = {x,T]), 

(10) \d:Hz)\<C^, \a\>2; 
(iii) there exists S > such that 

(11) \detdl^<!>{x,r])\>S. 
If we set 

.^2N ( y = V,<l>ix,7]) 

we can solve with respect to (x, ^) by the global inverse function theorem (see e.g. 
[21]), obtaining a mapping x defined by (x,^) = xiv^v)- As observed in [8j, x is a 
smooth bi-Lipschitz canonical transformation. This means that 

- X is a smooth diffeomorphism on R^"^; 

- both X and x~^ are Lipschitz continuous; 

- X preserves the symplectic form. That is, if (x, ^) = xiViV)) 

d d 

(13) 2_, dxj A d^j = 2^ dyj A dr]j. 

i=i i=i 

Observe that condition (TT3|) is equivalent to saying that the differential dxiy, rj), at 
every point {y, rf) G M^*^, as a linear map M?"^ — )■ M?'^, is represented by a symplectic 
matrix, i.e. belonging to the group 

Sp{d,R) := {A G GL{2d,R) : 'AJA = j} , 
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where 

. _ r /, 

In particular, when x is hnear, this means that x itself is represented by a symplectic 
matrix. 

We now address to the problem of the almost diagonalization of Fourier integral 
operators with a tame phase function. It turns out that, when the symbol is 
regular enough, the matrix representation of such an operator is almost diagonal 
with respect to a Gabor system. 

Theorem 11.11 in the introduction was the first precise result about the almost 
diagonalization of Fourier integral operators in time-frequency analysis: Assume 
that ^ is a tame phase function and g G 5(M'^). If a belongs to the Sobolev space 
W^{M.'^'^) (consisting of all distributions with essentially bounded derivatives up to 
order 2N), then there exists a constant Cn > such that 

(14) \{Trr{f,)g,rr{X)g)\<CM{x{t^)-X)-''', VA,/x eR'", 

where x is the canonical transformation generated by $. 

This statement proved in [51 Theorem 1]. For the extension of Theorem ll.il we 
next introduce the appropriate symbol class, li m & A^^^(R^'*), with Vg = Vs{x,r]), 
{x,r]) e ]R2d^ s G M, we consider the symbol class M^^(M.'^'^) of cr G ^'(IR^'^) 
satisfying 

lkllM°='^m2d^ = ||V*cr||r°o,i .^4ds = / sup \V^{z,C)\m{C)dC < oo. 

In time-frequency analysis this symbol class is just a special case of a modulation 
space [12] , in theory of pseudodifferential operators M^'^ is often referred to as a 
generalized Sjostrand class after 



For FIOs with symbols in M{^'^(R^'^) there is an almost diagonalization result 
similar to Theorem 11.11 

The window g will be chosen sufficiently regular, as follows. Given s > 0, let 

(15) iVGN, N>'- + d, and g e MI^^^^^R"), 

here v^n is a weight function defined on M'^. 

Theorem 3.2. Let m G Mv^iR^'^), s > 0, N, g satisfying flTSj) and such that 
Q{g,A) is a Parseval frame for L'^{R^). If the phase $ satisfies conditions (i), (ii) 
and a G M~'^(R2'i) then there exists H G W{Ll^){R?'^) such that 

(16) |(T7r(x,r/)(7,7r(x', V)^)| < H{r^' - VMx',v),x- V,<D(x',r/)), 
for every (x^rj), {x',r]') G M^^, with \\H\\w(lIJ < 1^11^^00.1 

Observe that Theorem 13.21 says we can control the coefficient matrix 
\{Tn{x,r])g,n{x',ri')g)\ by a function H of the difference (r/' — Vx^{x',ri),x — 
V^$(a;',?7)). It will be useful in the sequel to have a control of the matrix which 
depends on the difference (77' — X2{x,ri),x' — xi(x, r^)). This is achieved by adding 
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the assumption (iii) of |T] on the phase and by choosing another, even larger mod- 
ulation space as the symbol class. 

Theorem 3.3. Let s > 0, N, g satisfying ( fTSl) and assume that Q{g,A) is a 
Parseval frame for L'^{M.'^). If the phase $ is tame and a G M^^ {E?'^), then there 
exists C > such that 



(17) 



\(TA,),Mm\<c^^^., 



y\,jj, e 



o2(i 



Theorem [23] improves Theorem 1 1.1 1 in several respects: If s = 2N, then Wr 



Mr^^^JW^) m Prop. 6.7] and ^(M'^) C M^^^^^^ 



oo ('^^>2d^ 

2N\ 



C 



'^). Thus we obtain the same 
quality of off-diagonal decay for significantly larger classes of symbols and windows. 

The remainder of this section is devoted to set up the tools for the proof of these 
theorems. 

For z,w E M^"', let $2.2 the remainder in the second order Taylor expansion of 
the phase $, i.e., 



\a\=2 



<l>2,z{w) = 2 J2 / {l-t)d''<!>{z + tw)dt 



w 



a 



For a given window g, we set 
(18) ^Jw) 



^2ni^2,ziw)— 



g®g{w). 



The main technical work is to show that the set of windows ^ z possesses a joint 
time-frequency envelope. This property will allow us to replace the z-dependent 
family of windows '^ z by a single window in many estimates. 

Before proving the existence of a time-frequency envelope in Lemma 13.91 below, 
we first look at the phase factor e^'^**^'^ occurring in (fTSi) . 



Lemma 3.4. For every s eE, N efi, N > ^ + d, and "^ e S{R'^'^), we have 

(19) sup |V*e2-*^-|GCj,U.(^'')' 

with t'_4Ar and Vg being weight functions on R^"^. 
Proof. We proceed as in [8]. Let \I^ G 5(M^'^), then 



\V^e 



27rj$ 



2,z I 



U,W) 



Using the identity 

(1 - A^)^e~2.'C» = ^27rw)2^e-2"'^"', 

we integrate by parts and obtain 

1 



\V^e 



2"'*^-^(m,w) 



{2tiw) 



2N 



/ (1 - A^)^(e2'^'*2'^(«)T„^(C))e-2"'^-"'rfC 

JE2d 
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By means of Leibniz's formula the factor (1 — A^)^(e^'^**^'^''^''T„\l/(C)) can be ex- 
pressed further as 



^27ri$2,,(C) 



E Paid<^2AC))iTud^c^m, 



\a\ + \/5\<2N 

where Pa{d^2,ziC)) is a polynomial of derivatives of $2,2 of degree at most |a|. 

As a consequence of (fTO|l we have \pa{d^2,ziC))\ ^ C'o(C)^'"' foi' every z G M^'^ 
with a constant independent of z. Moreover, the assumption \E' G iS(]R^'^) yields 
that 

sup \Tud^.^\<CNAC-u)-' 

\I3\<2N 

for every i > 0. Consequently, 



\V^e 



2ni^ 



2,z 



U// /llj2d 



< 



< 



(27rw) 
1 






(C)^^^(C-n)-^rfC 



(27rw)2^ 
whenever i > AN + 2(i. Since —2N + s < — 2d by assumption, we obtain 

/ sup sup |Kte2"**2'^|(M,w)(u)-^^(u;)'ciw<oo, 

jR2d u€R2d 2GM2d 

whence (IT9l) follows. 



Q 



Remark 3.5. Notice that the two weights V-^n and Vg compensate each other in 
(11 9p . It is easy to see that some form of compensation is necessary. For example, if 
$(C) = ICP/2, then $2,2(0 — ICP/2 independent of z. A direct computation shows 
that the STFT o/e^''**^'^*^''^ with Gaussian window ^(C) = e"'^'''' has modulus, up to 
constants, e"'^!"^""^! /2^ ui,U2 G M^*^. This function belongs to L^'^^^„^(M^'^) if and 
only if —AN + s < —2d. This is in fact better than the condition —2N + s < —2d in 
the assumptions, due to the fact that here the derivatives of order > 1 of^2,z{0 (^''^^ 
bounded from above by {(), instead of by (C)^; O'S in the general case. This explains 
the presence of 2N instead of AN . In the general case the best upper bound for the 
derivatives c»/$2,z(C) ^■5 (C)^; so that the condition —2N + s < —2d should be sharp. 

Lemma 3.6. Let m G A^„(]R'^), v, i>,w be weight functions on M.'^, and {f^ : z G 
M^'^} C S'{R'^) be a set of distributions m S' . If 



(20) 

for given ip E S{ 
(21) 



sup \V^fz\ e L 



00,1 



^(R"'') 



), then for every h G M^^^J 



sup|r^2(M)|GLi,^„(M2'^). 

zm.'i 
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Proof. The proof of the multipher property is similar to the convolution property 
of modulation spaces in [6j. Since Vgf{x,ri) = (/ ■ Tj.g){ri) [T71 Lemma 3.1.1], we 
obtain the identity 

and consequently 

where the convolution is in the second variable rj. Now set F{x, rf) = sup^^^2d \V^fz{x, rf) \ 

sup^gK2d \{fz ■ T^<fy{r])\ and H{x,ri) = \V^h{x,r])\. Then 

sup \V^2{f,h){x,ri)\ < {F*r^H){x,r]). 



Finally, 

\\^*^^\\Ll^m= / / F{x,r]-C)H{x,C)dCw{x)m{r])dxdr] 

Jmd jRd j-^d 

— F{x,ri — Qiy{x)^^m{r] — QH{xX)uj{x)i'{x)v{Q dxdr]d( 



<^ II ni II trii 

^ \\r \\ ^oo.l Ln r 1 

In the last expression both norms are finite by assumption. D 

Corollary 3.7. Let s > 0, N e N, N > ^ + d and g e M^\^^^^^(M'^) and 

^ e >S(M2rf). Then the function "^.{w) = e^^**^,.^;^ (g, ^(^) defined as m (|T8l) 
satisfies 

sup \V^^^{u,,u,)\eL\^,X^"'). 



Proof. The symmetry property of the weight (x)"^^ (rj)'^^ implies that M^^^^^^^(M'^) 
is invariant under the Fourier transform, and thus g G M^^^^^^^(M'^). A tensor 
product argument then shows that g®g & M^(R^'^) with Vr(a;i,X2, 771,772) = 
(xi)^^ (0:2)^^ (771)^^ (7^2)'^. Since W^(xi,X2,7^i,7/2) > ((a:i,a;2))^^((77i,7;2))'^, we 
also obtain that g®ge Ml^^^^^i^^^"). 

By Lemma Sai we have sup^^jj^. \V^e''^'^^-^\ e C_'|^^„^(M4'^). 

We now apply Lemma [3.61 with /^ = e^'^**^^^ and h = g ® g, u = V4i\i,m = Vg, 
w = 1, and V = v^jy (observe that Vs is f4Ar-moderate, since s < AN). As a 
conclusion we obtain that sup2g]g2d jV^^E'^l C L|g,^^(R'^'^). D 

We recall the following pointwise inequality of the short-time Fourier transform 
[TTJ Lemma 11.3.3]. It is often useful when one needs to change window functions. 

Lemma 3.8. //^o,^i,7 e S{R'^) such that (7,^1) y^ and f E S'{R'^), then the 
inequality 

\VgJix,v)\ < jJ-y^i\V,J\ * |V^,o7l)(x,7/), 



holds pointwise for all {x, rj) G 



p2d 
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Lemma 3.9. Let s>0, N eN, N > ^ + d and g e M„\^^^^^(R"'). Then 
(22) sup \V^^^\ G WiLl^JiR'"). 



Proof. Let if E S(R'^^) such that \\^\\2 = 1- Using Lemma |3^ 

< \V^-^\ * ( sup \V<i,^ip\){Ui,U2). 



Since V^^H e 5(M^'^) C W{L\^J and \\snp,^^2,\V^M\\Ll^,^ < oo by Corol- 
lary [X3 the convolution relation Ll^^_^ * W{L\f^^J C W{L\^^J of [HI Theorem 
11.1.5] implies that 

II sup IK^^VPI WwiLl^^J < IKnWiLl^Jl ^nplV^M IIl1^„^ < oo, 
2GlR2d ^GK2d 

and so sup^g]g2d |V^*^^E'| G W{L\ 



®Vs 



fOO,l 



We next formulate an new characterization of the symbol classes M^^'^ and 
M^^ that is perfectly adapted to the investigation of Fourier integral operators 
with a tame phase $. 



Proposition 3.10. Lets>0,me M^^{R^'^), N, g be as in i^ and a e S'{M.'^'^). 

fOO,l / 



(i) Then the symbol a is in M^'^(M^'^), if and only if 



(23) II sup iV^^crlll^,""^! (184^) = / sup sup \Vxjj_,a{ui,U2)\m{u2) du2 < oo, 



with ||(t||^oo,i X II sup2p]R2d |V^^0"||| roo,i . 

In this case the function H{u2) = sup^^gj^d sup2g]g2d |V*^cr(ui,M2)| is inW{Ll^){R'^' 
(a) Likewise, a G M{^^(R^'^) if and only z/ sup2g]g2d |V*^cr| G L^^(]R'^'^) with 

(T MOO X SUp^cHjad Vij(,(T ro° 

II ll^"l®m II ^ZfcM I *z Ill-^l0m 

Proof. We detail the proof of case (i). Case (ii) is obtained similarly. Let "^ G 
e-^) with ll^lla = 1. 
Assume first that a G M^!^(R^'^). Then by Lemma [3.81 we have 

|V*^Cr(Mi,U2)| < \V^0-\ * \V^^'^\{ui,U2) < |V*0-| * SUp \V^^'^\{ui, U2). 



Set F{ui,U2) = sup2g]g2d |Vg,^\l'(ni,M2)|) so that 

\V^,cr{ui,U2)\ < {\V^cr\* F){ui,U2). 
Finally 

II sup I V*^cr| 11^00,1 = / sup sup \V^^a{ui,U2)\m{u2) du2 

zeR2d 1®™ JR2du^glRd^gK2d 

< II |V*o"| * F|| .00,1 

Ml I M Ij-. -- 



< II V^(T|| roo.i ||-F||ri < ||cr|Lroo,i ||-F||ri < OO, 
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where in the last step we used the independence of the weighted M°°'^-norni of the 
window |T71 Thm. 11.3.7] and Corollary 13.71 

To obtain the sharper estimate, set H{u2) = sup^_^^^dSup^^^2d\V<^^a{ui,U2)\, 
Fi{u2) = j^2dF{ui,U2)dui, and Gi{u2) = snp^^^^2d \Vq,a{ui,U2)\. Then the def- 
inition of Mi^'^(M2d) implies that Gi G ^^^(M^'i)^ and Lemma [33] implies that 
-Fi G W{Ll )(R^'^). With these definitions we obtain a pointwise estimate for H, 
namely 

H{w) = sup sup |V*^o"('Ui, w)| 



< / s\x\)\V^,a{ui,U2)\F{z — ui.w — U2)duidu2 

= {G,*F,){w). 

The convolution estimate (IHl) now yields that H < Gi * Fi G L]^ * W{L]j^) C 
ly (L^ J (M2"!), as claimed. 

Conversely, assume ( !23l) . Using Lemma 1X8) again, we deduce that 

|V;i;0r(Ml,M2)| < A y<i'M * \V^,'^z\{Ui,U2), 

and 

\{g (^7j){w)\^ dw = ||^1|2||?||2 = hWt 



is in fact a constant (depending on g). Using the involution f*{z) = f{—z), we 
continue with 

\Vii,a\{ui,U2) < \V^M * |Vvi,^2|(mi,M2) 

= \V^.(^\*\iV^.'^)*\iui,U2) 



< sup |V;i,^0-| * sup |(Vvi,^^)*| {Ui,U2) 

Since by Corollary 13.71 

sup liV^^^ri = sup \V^^^\ G Ll^^S^""), 



2GIR2'* 

we get 

||V^cr|Loo,i < II sup |V;^ ct||Loo,i II sup |V^^\E'|||ri < oo , 

and the proof is complete. D 

We now prove Theorems 13.21 and 13. 3[ 
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Proof of Theorem \3.2[ We know from [8, Equ. (39)] or [7, Equ. (3.3)] that 
I (T7r(x, r])g, 7i{x', r]')g)\ = \V^ ^^, ^^^^a {{x' , r]), (r/' - V^$(x', r]),x- V^$(a;', r] 
(24) < sup |Vv,,a((x',r/),(r/'-V,$(x',r/),x-V^<l'(x',r/)))|. 



We claim that the function that controls the off-diagonal decay of 
\{T7i{x,ri)g,'7i{x',r]')g)\ is exactly the function if (w) := sup„^g]j2d sup^gig2d |V^^cr(ui, w)| 
introduced in Proposition l3.10[ There we have already proved that H is in W{Ll^){R'^'^), 
and now ( !24l) implies that 

I {Tnix, r])g, 7r(x', 7]')g)\ < H{7]' - V^^x', 7]),x- V^$(x', t])). 

U 

A decay estimate in terms of the canonical transformation x can be obtained by 
imposing a stronger condition on the symbol class. 

Proof of Theorem \3.3[ We follow the pattern of the proof of Theorem 13. 2[ If $ is a 
tame phase function, then the argument of V^^a in (^^ can be estimated further 
as 

(25) \r]' -Vx^{x',r]) + \x -Vn^{x',r])\ > \r]' -X2{x,r])\ + \x' - Xi{x,r])\ 

by [U Lemma 3.1]. 

Set H{ui,U2) := sup2g]K2d |V*^cr(-Ui,'U2)|- By Proposition [3TT0l (ii). H E Lf^^^(R'^'^) 
and ||ii||roo < ||cr||M°o . Now, using (|2^ . 

\{T7i{x,r])g,7i{x',r]')g)\< sup H{ui,r]' -Vx^{x' ,t]),x -V^^{x' ,r])) 

< sup H{ui,r]' -Vx^{x',r]),x-Vr,^{x',ri)) 

MlgK2d 

[t]' - Vx$(x^ v),x- V-n^jx', 7])y 
{t]' - V^$(a;', 7]),x- V^$(x', t]))' 

< SUPh2 SUP^igR2d H{ui, U2)Vsiu2) 

~ {ri' - V^<l>(x', ri),x~ V^$(a;', ri))-' 

II Wll 

_ II "^10Vs 



< 



(V - Vx^{x', 7]),x- Vn^{x', 7])y 



where in the last inequality we have used f l25|) . D 

4. GaBOR MULTIPLIERS 

We now introduce a type of Gabor multipliers that are tailored to the canonical 
transformation x of a FIG. In the case of the identity map x = id we get the 
standard Gabor multipliers studied, e.g., in [T^ [20]. 

Since in general the map x does not preserve a given lattice A, we first replace 
X by a mapping x' that enjoys x'(A) C A. 
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We define the integer part of a vector y = {yi, ..., y2d) as 

Vy\ = (Li/iJ>---> Vy2d\), 

i.e., by taking the integer parts of its components. 

Now consider a lattice of the form A = AT?''- with A G GL{2d, M). Given a map 
X '■ K^'^ -^ IR^°', we approximate x by a map x' : A — > A defined by 

x'{\) = A[A-'xi\)\. 
Observe that \y — [y\\ < \/2d and therefore 

|x(A)-x'(A)|<y2dP||. 

The almost diagonalization of Theorem 13.31 can now be formulated in terms of x' 
as follows: using the inequality (x) < v^(x + y)(i/) we have 

(26) |(T7r(^)(7,7r(A)(7)| < C^Ti^^d\\A\f^ {^{^i) - A)-2^ VA,/i G A. 
We collect here some properties which will be used in the sequel. 

Lemma 4.1. Let v{z) = Vs{z) = {zY , z G M^'^, s > 0, and m G A^^(]R^'^). Then 

(27) vox^v, vox' ^v, 
and 

(28) moxGA^„(M^'^), m o x' G A1^(R^^). 
Proof. Since x is Lipschitz continuous, we have 

\x{z)\ < |X(0)| + Ix(^) - X(0)| < |x(0)| + C\z\ < C'{1 + \z\) Wz G R'' . 

By applying the same argument to x~^) which is also Lipschitz continuous, we 
obtain 1 + \x{z)\ x 1 + \z\. Since we also have x'i^) = x{^) + ^^(1); both the 
estimates in f l271) follow. 

Concerning fl25]) . observe that 

m{xiz+C)) = m{xiz+C)-xiz)+xiz)) < v{xiz+C)-xiz))Tn{xiz)) < v{C)m{xiz)), 

where in the last step we used the fact that \x{z + () — x(^)l ^ ICI- This proves 
the first formula in (128|) . The second formula is obtained similarly, because x'{^) = 
X{z) + C(l), so that x'{z + C) - x'{z) = x{z + C) - x{z) + C(l), and we still have 

i + lx'(^ + C)-x'(^)l<i + IC|. D 

Remark 4.2. Notice that the above proposition holds only for weights with polyno- 
mial growth, but not for weights with super-polynomial growth, e.g., v{z) = e"^'^' , 
0<6<l,a>0. This is related to the fact that the weights Vg satisfy the doubling 
condition Vs{2z) x Vs{z), whereas weights with super-exponential growth do not. 

Now, let Q{g, A) be a Gabor system, and x? x' be the canonical transformations 
of an FIO. Given a sequence a = {0'u),y^A' ^^ define (formally) the Gabor multiplier 

MJ = M^'^'^f = J2<^x {f,7T{X)g)7T{x'{X))9. 
AeA 
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In order to give a precise meaning to this definition, we need to study the conver- 
gence of the above series. 

Lemma 4.3. Let g E M^{R'^), m E A^^(R2^), andl <p < oo. If h = {hx)xeA G 
Co,. (A), then 



(29) \\J2hxHxW)9\\M^^<c[j2\^xfm{x'i\)y 



AeA agA 



The series YlxeA^>^'^ix'W)9 converges unconditionally in M^(R ) forp < oo and 
weak* unconditionally for p = oo. 

Proof. We assume p < oo and leave the modification for the case p = oo to the 
reader. 

Since the finite sequences are norm-dense in C(^) ^^ p < oo (and weak*-dense 
for p = oo), it suffices to show (!29|) for finite sequences. For a finite set F C A we 
set F = 'x'{F). Since g E Mj(]R'^), the synthesis operator {cq,} h- X]aeA'^Q^('^)fl' i^ 
bounded from ^(A) to M^ (M'^) [HI Thm. 12.2.4], so that 

II ^/iA7r(x'(A))c/||M^ = II 5Z ( 5Z ^A)7r(a)^||M^ 

(30) <(E| H h> ^ 

aeF Aex' ^(a) 



mfa)^ 



Since x ^ is Lipschitz continuous, x' is "uniformly almost injective", in the sense 
that supogy^^ #x'^ ({<^}) < oo. Hence the last expression in (l30!l is bounded by 



P 



^(E E i^aIm«)^) ' = (Ei^^imx'(a)) 

By density (12^ holds for all (/ia) G Cov'(^)- The unconditional convergence 
follows from the norm estimate. D 

Proposition 4.4. Let g E M^2{M.'^), m,m E Mv(R'^'^), and 1 < p < oo. If a = 

{ci\)\eA is a sequence in £^(A), then the Gabor multiplier M^ = M^ '^'^ is bounded 
from M^ , {R'^) to M^{R'^). Its operator norm is bounded by ||Ma||jv/p ^m'' ^ 



Proof Since ^ E M^2{R'^) and g E M^^iR'^), the coefficient operator / ^ 
(/, 7r(A)^) is bounded from M^(M'^) to l^J^iA), see [HI Thm. 12.2.3.]. Hence, if 

m m 

a = (aA)AGA e ^^(A) and / E M^(R'^), then the sequence (aA(/, 7r(A)5f))AeA be- 

longs to Cox'(A)- Therefore, by LemmaiJlthe series X]aga'^a(/, vi"(A)5f)7r(x'(A))fi' 
converges unconditionally in M^(]R'^). The estimate for the operator norm follows 
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from 

aga AeA 

<||*a||..(El(/.^(%>K^)T" 

AeA ^ ' 

^ ll«l|£^ll/llAf^ ,• 



Corollary 4.5. //a G £^(A) for s e^ and g G 5(M'^), then M^''^'^ maps 5(M'^) 
mto>S(M'^). 

Proof. Since w^v o x' x w^r by Lemma SH M^^^^,/^^(]R'^) = Mp^_^(M'^). By Propo- 
sition HaiM^ maps MP^_^{R^) to M^{R^). Since 5(M'^) = n7v>o ^^°^ (^'^) ^Y ©, 
Ma maps 5(M"') = n7v>o ^^°^-. (^'') i^^o S{R'^). ~ D 



5. Approximation of FIOs by Gabor Multipliers 

Our next aim is to approximate a FIO T with a tame phase by Gabor multipliers 
associated to a Parseval frame Q{g,A). First we will derive a representation 



(31) T = J2^iu)M,^, 

with convergence in several operator norms. To find the candidate symbols a^, we 
argue as in [20] . 

Using the commutation relations M^T^ = e'^^'^^^T^M^j, we can write 

(32) 7i"(x'(/^) + 1^) = c^,^7r(z/)7r(x'(/i)), for z/, /i G A , 
with \c,^,i_i\ = 1. We will show that the choice of 

(33) a^(/i) = c^^f,{Tn{iJ,)g, 7r(x'(/i) + iy)g) 

leads to the formal representation (13T|) . 

Since ^((7, A) is a Parseval frame, / and Tf possess the Gabor expansions 

/ = 5^(/, rr{f,)g)iT{f,)g, and T/ = ^^(T/, 7r(A)^)7r(A)(7. 
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Setting z/ = A — x'(/i) G A, we can write 



fj, X 

= X^ X^(/, 7r(/i)^) (Tn{fi)g, 7r(x'(/i) + u)gjn{x'{f^) + i^)g 

/I V 

= X] X^ CvM^ 7r(/i)^) ^T7r(/i)^, 7r(x'(^) + v)g^ 7r(z/)7r(x'(/i))^ 
= X^ 7r(z/) ^ Cj.,^(/, 7r(/i)^) {T7r{fi)g, 7r(x (/x) + v)g)Tx{x{^i))g 
(34) =5^7r(z/)M.J. 



The following result gives a precise meaning to the above computation for test 
functions in the Schwartz class i5(M'^). 

Proposition 5.1. Let Q(g,A) be a Parseval frame and g E S{W^) \ {0}. 

(i) If |ai^(yu)| < C(l + |/i| + |z/|)^/or some N > 0, then the series Xl;/eA^('^)^a^ 
converges unconditionally in the strong topology of C{S{M.'^),S'{M.'^)) and defines 
a continuous operator from 5(R°') to iS'(]R^) (S'lMf^) is always endowed with the 
weak* topology). 

(a) Let A be a continuous operator from 5(M°') to 5'(M°'), and choose ay(/i), I'yfJ.E 
A as 

a^(/i) = c^,^(A7r(/i)5f,7r(x'(Ai) + iy)g) 

with the constant Cy^^ as in fl52]) . Then the grow estimates 

(35) |a,(/i)|<C(l + |/i| + |i/|)^, 

are satisfied for some constants C,N > depending on A and g. 

Furthermore, A can be represented as the sum '^i,^j^'n'{i')Ma^ of shifted Gabor 
multipliers in the strong operator topology of C{S{W^),S'{W^)). 

Proof, (i) By Corollary 14.51 the Gabor multiplier 

AteA 

maps 5(R'^) to 5(R'^). We need to show that the series J^^eA T^i.^)M^J, f ^ '5(M^), 
converges unconditionally in 5'(R'^). 
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Let h G iS(M'^) and choose s and / large enough, precisely s > N + 2d and 

l> N + s + 2d. Then 



< 



5^ 5^ |a.(/u) 1 1 (/, 7r(/i)(7) 1 1 (7r(z/ + x'(/x))(7, /i) I 



1/ fi 



(36) ^ E E(l + l/^l + IH)"' (1 + l/^l))"'(l + IX (^)l)^(l + IH)-^ < oo. 



1/ /i 



In the last step we used the assumption |a,^(/i)| < (7(1 + |/i| + |i^|) and the fact that 
1 + |x'(/i)| X 1 + |/i| from Lemma [4.11 Then the double series converges because of 
our choice s > N + 2d and I > N + s + 2d. 

(ii) The proof is similar to that of Proposition 5 (ii) of [20] • Let K G 5'(]R^'^) 
be the Schwartz kernel of the operator A. Writing fi = {fii,fi2), ^ = {^i, ^2) ^ ^'^'^ 
and x'(/^) = ixM^xM) ^ ^*'^, we have 

\aM\ = \{ATT{fi)g,TT{x'{fJ') + Jy)9)\ 



= \{K,n{x'{fi) + iy)g^7f{fi)g)\ 

= \Vg(x,gK (x'lifi) + iyi,fii,xM + ^2, -^2)1 

<c{l + W{^^)\ + \^^\ + \v\r 

for some C, A^ > 0. In the last step we used the fact that the STFT of a tempered 
distribution grows at most polynomially. Since 1 + |x'(yu)| ^ 1 + |yu| "we see that the 
sequence Sii,{fi) satisfies the polynomial growth condition ( l35l) . It follows then from 
part (i) that the series J2,yeA'^i^)^a^ converges in £(5(R°'),5'(]R'^)). Furthermore, 
its sum must coincide with A, because the formal computation flM|) was justified 
in part (i) and therefore 

{J2 A^)M.J, h) = {Af, h), V/, h G SiR'). 

u 

(One needs to interchange two summations, which is possible by fl36l) ). D 

Remark 5.2. Proposition 15.11 extends to more general windows provided that the 
sequences a^, satisfy stronger estimates. 

For example, if sup^^^g^^ |aj^(/i)| = C < 00, then one may use windows g G 
M\R'^). In fact, for /, A G M^W^), (EE]) can be modified to yield 

<C5^|(/,vr(/.)^)|5^|(7r(z. + x'(/i)k,/^)l 

fj. V 

<C'\\n\MA\h\\.n <00. 



where in the last step we have used the characterization ([5]) of M^ (M'^) (with m = 1 
and p = 1). 
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We now consider the convergence of series of Gabor multipliers on modulation 
spaces. In what follows the space A^J^(M'^) denotes the closure of the Schwartz 
class with respect to the M^-norm. Whereas Gabor expansions converge only 
weak* on M^{W'-), they are norm convergent of Ai'^iW^). 

Proposition 5.3. LetQ{g,A) be a Parseval frame and g G M^(]R'^). If the sequence 
of symbols a^ satisfies 

y ||a,^||£oot>(z/) < oo , 

then the series J2u '^{^)^a^ converges in C{M^ ,, M^) for every p G [1, oo) and v- 
moderate weight m. Ifp = oo, the series J^u 7r(z^)Ma^ converges in C{J^'^^^, , Ai'^) . 

Proof. By Proposition l4.4l each operator Mg,^ is bounded from M^ , (M.'^) to M^(M.^) 
with the norm being dominated by ||ai,||^oo. Furthermore, \\'n'{i^)\\M^^M^ ^ '"{^)- 
Hence 



J2 IK('^)^aJU/^^^^,^Af^ < Yl ll^(^)IU<^A/fJl^aJ|Af^^^,^Af^ 
u V 

< ^ ||a,,||^oot;(z/) < oo. 

V 

This gives the desired conclusion. D 

We return to the study of Fourier integral operators. The following approxima- 
tion theorem is the main result of our work. 

Theorem 5.4. Let Q{g, A) be a Parseval frame and g G iS(R'^). Let T be a Fourier 
integral operator with a tame phase and with the associated lattice transformation 
x' , and define the multiplier symbol as 

^M = c^,f,{TTT{iJ,)g, 7r(x'(/i) + i^)g) ■ 

(i) If a e W^{^^'^) for N > d and if < r < 2N - 2d, then the series 
J2u&A^i^)M^. converges to T in C{M^^^,,Mp^) and in C{M'^„^, , M^) , for every 
p G [1, oo), and Vr-moderate weight m (vr{z) = {zY , z G M^'^j, and 



\u\<L 



a^llAf' ,-^MFn 



mox 



(ii) If the symbol a is m M^^^{^^'^) with s > 2d, then T = Ez.eA ^('^)^a. with 
the error estimate 



\T- J2^i^)M.jMr' ,^Mr^<U-^"'- 

* ^ mov' 



mox 

\v\<L 



Proof, (i) We first estimate the magnitude of the multiplier symbols a^,: Theo- 
rem [TTTl written with x' as in fl2B]) . implies that 

|a.(/.)| = |(r7r(/x)^,vr(x'(/i) + v)g)\ < (x'(/i) - (x'(/^) + ^))-'^ = (^)~'^, 
for every yU G A. Hence 

lla^lloo < (^)-^^ 
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and 

By Proposition 15.31 the series Xlj^ '^('^)^a^ converges in the operator norm from 
M^ ,{R'^) — >• M^(R'^) for every 1 < p < oo. Since the series converges to T in 
£(5(]R'^),5'(M'^)) by Proposition 15. 11 the sum must be identical to T. 
For the error estimate we observe that 

W\<L W\>L 

\u\>L 
^ V^ I \r-2N ^ Tr+2d-2N 

\u\>L 

(ii) is proved similarly by using the decay estimate of Theorem 13.31 instead of 
Theorem 11.11 D 

As a byproduct of Theorem 15.41 we obtain an alternative proof of ^ Theorem 

4.1]: 



Corollary 5.5. Under the assumptions of Theorem \5.4l the Fourier integral oper- 
ator T is a bounded operator from M^ ,{R'^) to M^{M.'^), simultaneously for every 
l<p<oo, and from M^^^,{R'^) to M^{R'^). 

Remark 5.6. So far we have used without loss of generality that Q{g,A) is a 
tight Gabor frame. IfQ{g,A) is an arbitrary frame with g G M^{M.'^), then there 
exists a dual window 7 G M^lW^) such that every f G M^(]R'^), for 1 < p < 00 
and m G A^t,(M^'^), possesses the Gabor expansion f = X]aga(/' ^('^)5')^('^)7 '"^^^^ 
convergence in Mp^IW^). For the existence of a dual window in M^{R^) see [2T], 
for the result on Gabor expansions see [17J. All results of Sections 2-5 carry over 
to general Gabor frames by polarization. 

6. Further Remarks 

We conclude with an example and some remarks on how Gabor frames are trans- 
formed under FIO. 

We first compute explicitly the multiplier symbol ( 155|1 for the dilation operator 
in dimension d = 1 

Dsfix) = f{sx) = [ e2--V'(r/) dr] s>0. 

This is a Fourier integral operator with symbol a = 1 and phase $(x,?7) = sxrj. 
Using (IT2|) . the canonical transformation x is calculated to be 

x{y,v) = {y/s,sri). 



20 ELENA CORDERO, KARLHEINZ GROCHENIG AND FABIO NICOLA 

For the lattice A = aZ x /3Z with a/3 < 1 and fi = {ak, (3i), k,i & Z, we then have 
x'ifJ-) = {a[k/s\, f3[si\). Consequently, the symbols are 

aM = (^7r(/i)^,7r(x'(/i) + '^)9) 

^ g2™/3LtJ£' f ^2.^pisi-lsH-ntg^,^ _ ^k)-g{t -a[-\~ ak') dt. 
Jr -s 

Now, we consider the case of the window function g(t) = e~'^* . A straightforward, 
but lengthy computation (as in fiO\) gives 

Vs^ + 1 

The representation of the dilation operator D^ by a (sum of) modified Gabor 
multipliers might be of interest for the numerical approximation of Dg within the 
context of time-frequency analysis, when one is forced to use Gabor frames. Similar 
formulas can be worked out for general metaplectic operators as well. 

More on Gabor frames: Assume that Q{g,A) is a frame for L^(M'^) and that T 
is a FIO with canonical transformation x satisfying the standard conditions. Then 
the transformed system TQ{g, A) = {T7i{X)g : A G A} is a frame, if and only if T is 
invertible. However, if we fix the window g and only warp the time-frequency space 
with X, then we obtain the set Q{g, x(^)) = {''^{xW)g '■ -^ ^ ^}- It is an interesting 
problem to determine when the transformed Gabor system is still a frame. 

If g{t) = e"'^* in dimension d = 1, then Q{g,x{-A)) is a frame, if and only if 
the lower Beurling density d~{x{A.)) > 1. This follows from the characterization of 
(non-uniform) Gabor frames by Lyubarskii and Seip [25l l26] . 

On the other hand, if T = D^ is the dilation operator, A = aZ x /3Z, and g is 
a window with compact support, suppt/ C [—A, A], say, then x(^) = ~^ ^ sf3Z 
is again a lattice. Choosing - > 2 A, then different translates g{t — s^^ak), A; G Z 
have disjoint support, and Q{g, x(A)) cannot be a frame. Thus the frame property 
of Q{g, x(^)) is subtle. 

Based on coorbit theory |13] one can formulate a qualitative result. Recall that a 
(non-uniform) set A C M."^^ is relatively separated if max^j^^^^ cardAn(/c + [0, l]^'^) < 
oo, and A is called 5-dense for 5 > 0, if [J^^jyBs{X) = M.'^'^. 

The results in [13] assert that for every g G M^ (W^) there exists aS > depending 
only on g, such that every relatively separated (5-dense set A generates a frame 
g{g,A). We fix c/ G M\R'^) and 6 = 6{g) > 0. 

Now let T be an FIO with canonical transformation x with Lipschitz constant 
L. //AC M^'^ is 6/L-dense, then Q{g,x{-A)) is a frame. 

To see this, observe that 

inf|z-x(A)| = inf|x(x-'W)-x(A)| 
agA AeA 

< L ini \x'\z) - X\ < L- = 5. 
~ AeA L 

Consequently x(A) is 5-dense, and so Q{g,x{-^)) is a frame. 
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